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O ! ABSTRACT 

<N 

— h " Context. The mass of the neutron star in Vela X-1 has been found to be more massive than the canonical 1.5 Mq. This result relies 

on the assumption that the amplitude of the optical component's measured radial velocity curve is not seriously affected by the 
interactions in the system. 

Aims. Our aim is to explore the effect on the radial velocity curve caused by surface motions excited by tidal interactions. 
Methods. We use a calculation from first principles that involves solving the equations of motion of a Lagrangian grid of surface 
elements. The velocities on the visible surface of the star are projected along the line-of-sight to the observer to obtain the Doppler 
shifts which are applied to the local line-profiles, which are then combined to obtain the absorption-line profile in the observer's 
reference frame. The centroid of the line-profiles for different orbital phases is then measured and a simulated RV curve constructed. 
Models are run for the "standard" (v sini=\\6 km/s) and "slow" (56 km/s) supergiant rotation velocities. 

Results. The surface velocity field is complex and includes fast, small-spatial scale structures. It leads to strong variability in the 
photospheric line profiles which, in turn, causes significant deviations from a Keplerian RV curve. The peak-to-peak amplitudes 
of model RV curves are in all cases larger than the amplitude of the orbital motion. Keplerian fits to RV curves obtained with the 
"standard" rotation velocity imply m„, >1.7 M Q . However, a similar analysis of the "slow" rotational velocity models allows for 
m„ s ~1.5 M Q . Thus, the stellar rotation plays an important role in determining the characteristics of the perturbed RV curve. 
Conclusions. Given the observational uncertainty in GP Vel's projected rotation velocity and the strong perturbations seen in the 
published and the model RV curves, we are unable to rule out a small (~1.5 M Q ) mass for the neutron star companion. 
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' 1. Introduction processes (Woosley & Heger 2007; Timmes et al. 1996). Van 

_l ' , ... „ ,. w , „ . ,.„„„„. .„ den Heuvel (2004) has argued that the GP Vel system provides 

. The high mass X-ray binary system Vela X- 1 (4U 0900-40) con- direct evidence ^ there ^ ^ least Qne of neutnm ^ 

^ . sists of a pulsar (P, =283 s) in an eccentric (e=0 0^ P 6 =8 96 fc indeed fe i& such , masfj 

. ■ d) orbit around the B0.5-supergiant star HD 77581 =GP Vel. _ . , . . , . ,. 

_j ■ t, .. , • . . • ... . • , j , a ■ a The procedure lor determining the neutton star mass in bi- 

^ Particular interest in this system is sparked by the mass derived F , , , , ,. , , . ,™ 

O c . . \* i o nx • Ji- ■ ■ c .1 nary systems such as GP Vel is based on the radial velocity (RV) 

lor the neutron star, M„, ~1.8 M , since this is significantly J J ..... . , r , ,. , . .. . 

PnJ , • i 1 c A/r ,. .... . au curve which is derived from photospheric absorption lines ob- 

v N larger than the canonical 1 .5 M upper limit that is predicted by ,., r , • , • t ,• , ■ ■ • 

t— I ■ . j j . . .• / . . „ , ,. served in the spectrum of the optical companion. Implied in this 

. . , the standard neutron star equation of state as well as that which is , ■ , • , , ™ , • • , 

>j-., f „ j , e 1f , . approach is the assumption that the RV variations truly repre- 

predicted from supernova models for the newly formed neutron , ... . n , T , r ^^i, , , 

' ^ ' stars sent orbital motion of the star. In the case of GP Vel, the 

T , , , f . . , , validity of this assumption is highly questionable because it has 

rN . The maximum possible mass of a neutton star depends on the , , J , , F . , , . ^ , r , 

5— i .• e <- ^ f .i, u a /t **• long been known that the photospheric absorptions of GP Vel 

3 ■ equation ol state ol the ultradense matter in its interior (Lattimer f . ,. K , .,. ,, , ,. 

& Prakash 2007; Page & Reddy 2006). For nuclear densities ™ de , r g° P r ™ ent llne P ro u file variability. Specifically, the lines 

t„ia14 -3 4.x. ,.• e t t ■ " f." a t\. develop asymmetries over the orbital cycle which lead to system- 

~3xl0 g cm , the equation ol state is a soft one, and the . , F . { . , , „, , , J . r „ , . 

■ , C w T-, , i ... .1 .. f atic deviations of the measured RV data points from a Keplerian 

maximum mass is ~ 1.5 M . For larger densities, the equation of T . . ,. , , . ,. F 

. . • „ , ■ j * , . j .. e RV curve. It is believed that the line-profile variability is associ- 

state is a stiff one, and depending on the adopted equation of ,. , , • , • 

.... ,. .. , ■ o n i\/r iv i ated with tidal forcing, non-radial pulsations and other mterac- 

state, the upper mass limit may be as high as 2.9 M (Kalogera . _ . , , . fe ' , , 

„ D inn^ t-u • ■* e * * ■ u- * tion effects in the binary system (van Paradiis et al. 1977b; van 

& Baym 1996). The majority of neutron stars in binary systems Tjr . ... . _/ < , „ . „ , 

do indeed have m ns <1.5 M (Thorsett & Chakrabarty 1999; ^erkwijk et al. 1995; Barziv et al. 2001, Quaintrell et al. 2003). 

Schwab et al. 2010; Lattimer & Prakash 2010). On the other Thls ralse , s * e question of whether these systematic deviations 

• . e 01V . . . • e ■ , may not lead to an artificially large RV curve amplitude and 

hand, the existence of m„ s ~2 M neutton stars is now fairly . J . <■ 

ii * i,t i. a m. * . i oaian nu u v- t. i hence, an overestimate of the neutron star mass, 

well established (Demorest et al. 2010), although it is not clear 
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whether such high-mass neutron stars are created at birth in the The tidal effects are clearl y present, as shown by the opti- 

supernova event or are the consequence of subsequent accretion cal h g ht curve whlch displays ellipsoidal variations with a full 

amplitude ~0. 1 mag indicating that the star is strongly distorted 

Send offprint requests to: G. Koenigsberger by the neutron star's gravitational field (Jones & Liller 1973; 
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Zuiderwijk et al. 1977; Tjemkes et al. 1986). Van Paradijs et al. 
(1977b) explored the effects of the deformation of the star on the 
RV curves. They assumed a circular orbit and synchronous ro- 
tation and computed the shape of photospheric absorption lines 
assuming a non-uniform temperature distribution. However, the 
rotation of GP Vel is clearly slower than synchronous, and it 
cannot be assumed to have a simple Roche geometry. Indeed, 
Tjemkes et al. (1986) showed that a model for the ellipsoidal 
variations based on the assumption of an "equilibrium tide" does 
not adequately reproduce the shape of the observed light curve. 
They suggested that non-linear effects might be important, but 
full non-linear calculations for tidally interacting stars are only 
now becoming available (Weinberg et al. 201 1). 

We have developed a 2D calculation, the TIDESQ code, that 
provides the time-dependent shape of the stellar surface and its 
surface velocity field for the general case of an elliptic orbit and 
asynchronous rotation. Using the derived velocity field, the line- 
profile variability is computed (Moreno & Koenigsberger 1999; 
Moreno, Koenigsberger & Toledano 2005). The method, though 
limited to the analysis of the surface layer, provides insight into 
the non-linear effects that appear in binary stars that are not in 
an equilibrium configuration. In particular, it allows the analysis 
of highly eccentric and asynchronous systems. Here, tidal flows 
are present in forcing regimes that are far from equilibrium con- 
figurations, and linear approximations are inapplicable. In this 
paper we use this model to explore the effects on the absorption 
line-profiles produced by the tidal flows on the B-supergiant's 
surface, and the resulting deformation of the RV curve. 

In Section of 2 of this paper we summarize the method for 
calculating the tidally-perturbed RV curves; Section 3 describes 
the method for chosing the stellar and binary parameters; Section 
4 contains the results; and Section 5 lists the conclusions. 



2. Calculation of tidally-perturbed RV curves 

2.1. The TIDES code calculations 

Our method consists of computing the motion of a Lagrangian 
grid of surface elements distributed along a series of co-latitudes 
covering the surface of the star with mass m\ as it is perturbed 
by its companion of mass ;«2, which is assumed to be a point 
source. The main stellar body below the perturbed layer is as- 
sumed to have uniform rotation. The equations of motion that 
are solved for the set of surface elements include the gravita- 
tional fields of m\ and nt2, the Coriolis force, the centrifugal 
force, and gas pressure. The motions of all surface elements are 
coupled through the viscous stresses included in the equations 
of motion. The surface layer is coupled to the interior body of 
the star also through the viscosity. The simultaneous solution of 
the equations of motion for all surface elements yields values of 
the radial and azimuthal velocity fields over the stellar surface, 
v r and v^, respectively. 

The calculated velocity field is then projected along the line- 
of-sight to the observer to obtain the Doppler shifts required 
to produce the integrated photospheric absorption line profiles. 
We assume a Gaussian shape for the local profiles^ and a limb- 
darkening law of the form s(6) = (1 — u + ucosO), with m=0.6 
as in the Milne-Eddington approximation. Full details of the 
model are given in Moreno & Koenigsberger (1999), Toledano 



1 Tidal Interactions with Dissipation of Energy through Shear 

2 A discussion on the effect on the results of the intrinsic line-profile 
shape is given in Harrington et al. (2009); the use of Gaussians is sup- 
ported by the results of Landstreet et al. (2009). 



et al. (2007) and Moreno et al. (2005, 2011). An application of 
the model to the B-type binary a Vir (Spica) may be found in 
Harrington et al. (2009). 

It is worth noting that the TIDES model computes the sur- 
face motions of the optical component from first principles and 
allowing for non-linear effects. However, there are currently two 
simplifications in the method. The first simplification is that the 
equations of motion are solved only for a thin surface layer, in- 
stead of for the entire star, which has the disadvantage that the 
perturbations from deeper layers in the star are neglected. This 
includes the excitation of non-radial pulsation modes which have 
been shown to cause oscillations in the RV curve (Willems & 
Aerts 2002). On the other hand, its response is representative of 
the effects that are produced on the outer stellar layer, which is 
the one that is most strongly affected by the tidal forces (see, 
for example, Dolginov & Smel'chakova 1992). The second sim- 
plification is that motions in the polar direction are suppressed, 
allowing only motions in the radial and azimuthal directions. In 
addition, we neglect the effects of possible temperature varia- 
tions across the stellar surface. 

The benefits of the model are: 1) we make no a priori as- 
sumption regarding the mathematical formulation of the tidal 
flow structure since we derive the velocity field from first prin- 
ciples; 2) the method is not limited to slow stellar rotation rates 
nor to small orbital eccentricities; and 3) it is computationally 
inexpensive. 

The parameters that are needed for the calculation are: the 
orbital period, P, eccentricity, e, inclination, i, and argument of 
periastron, <±> per \ the stellar masses, m\ and ma, the radius of the 
primary star, R\, its equatorial rotation speed, v rot , the polytropic 
indejfl n, and kinematical viscosity of its material, v. In addition, 
the code requires the depth of the surface layer, dR/R\ and the 
number surface elements for which the equations of motion are 
to be solved. The latter is specified by the number of longitudes 
into which the equator is divided and the number of latitudes 
between the equator and the polar region. We have done a thor- 
ough investigation of the TIDES code behavior (Harrington et 
al., in preparation) and find that 500 partitions in the azimuthal 
direction at the equator and 20 latitudes is sufficient to resolve 
the small-scale structure. This implies approximately 6800 sur- 
face elements^ in the semi-hemisphere above and including the 
equator. Since the axis of stellar rotation is assumed to be paral- 
lel to that of the orbital motion, the perturbations on the northern 
and southern hemispheres are assumed to be symmetric. Table 
1 gives a description of the input parameters and lists the values 
for those parameters which were held constant throughout all the 
calculations in this paper. 

The line-profile calculation is performed after the initial tran- 
sitory phase of the calculation has damped down. In the case of 
GP Vel, the steady state is attained after ~20 orbital cycles, and 
the line-profile calculation is performed at 30 orbital cycles. 

2.2. RV measurements 

For each given set of input parameters, the output of the TIDES 
code yields two line-profile files, one containing the line-profiles 
arising in the tidally perturbed surface and the second containing 
the profiles arising from the unperturbed, rigidly-rotating sur- 
face. The lines are measured to obtain their centroid. The cen- 



3 The polytropic index of a radiative supergiant star is n >3, 
Schwarzschild, 1958, p.258. 

4 The number of azimuthal partitions decreases from the equatorial 
latitude to the polar region. 
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Table 1. Description of input parameters for the TIDES code. 2.4. A note on orbital phases 



Parameter 


Description 


Values used 


Mi (Mo) 


mass of star 1 


see Table 3 


M 2 (Mo) 


mass of star 2 


see Table 3 


Ri (R G ) 


stellar radius at equilibrium 


see Table 3 


P orb (days) 


orbital period 


8.964368 d 


e 


orbital eccentricity 


0.0898 


i (deg) 


inclination of orbital plane 


78.8, 85.9 


u> per (deg) 


argument of periastron 


332 


A> 


asynchronicity parameter 


see Table 3 


v (R 2 G /day) 


kinematical viscosity 


0.12, 0.22 


a 


absorption-line depth 


0.7 


n 


polytropic index 


3.0 


k 


line-profile broadening param 


30 




num. segments along equator 


500 


N to 


number of latitudes 


20 



troid was computed through the weighted sum over all wave- 
lengths, At, as follows: 



A ce 



£/t.-(C,- - Id V2 
KG-/,-) 3 / 2 



(1) 



where the summation is carried out between the two limits 
within which the absorption line lies. In practice, the limits are 
the points where the absorption meets the continuum level. /, and 
C, are the line and continuum intensities, respectively. This defi- 
nition of the centroid is the same as that used by the 'e' function 
in the IRAI0 subroutine splot . 

The RVs were also measured through a Gaussian fit to the 
line profiles, using the IDL routine GAUSSFIT. Although a 
Gaussian fit is a very poor approximation to the actual shape 
of the lines, the derived RVs were in general very similar to 
those obtained with the flux-weighted centroid method. In the 
remainder of this paper we use the flux- weighted centroids. 

2.3. Semi-amplitude of the RV curve 

The radial velocities that were obtained from the model line pro- 
files were used to construct the RV curves. These curves were fit 
with the function that characterizes the orbital motion; i.e., the 
function describing the Keplerian orbifj 



V r — Vo + K\ [e COS U) per + COs(0 + <x)p er )\ 

where 



Ki = 



2n a\ sin i 
~P (1 -e 2 ) 1 ' 2 



(2) 



(3) 



is the semi-amplitude of the RV curve, a\ is the semi-major axis 
of mi's orbit, Vo is a constant vertical offset (the "gamma" ve- 
locity), and 6 is the true anomaly. 

The fit was performed using the Powell minimization method 
in an IDL script]. The free parameters for the fit were a \ sin i and 
Vq. The fitted value of K\ was then computed using eq. (3). 



5 Image Reduction and Analysis Facility, distributed by NOAO 

6 Binnendijk, 1960,p. 149-151 

7 The POWELL algorithm as described in Press et al. 1992, Section 
10.5 



In order to be able to compare the theoretical models with the 
observational data, the conventions for setting orbital phase 0=0 
need to be specified. In this paper, we define 0=0 at periastron 
passage. In most observational investigations <t> b s =Q is defined 
at the orbital latitude 90° from the line of nodes. This definition 
is such that <p bs=Q at the midpoint of the X-ray eclipse (Bildsten 
et al. 1997, van Kerkwijk et al. 1995, Barziv et al. 2001), and 
periastron passage occurs at O /, S =O.17. Hence, O /, S =0-O.17. 



3. Observationally constrained parameters of GP 
Vel 

The spectral type and class of the optical component in GP Vel 
is given by Morgan et al. (1955) as B0.5Ib. Howarth et al. (1997) 
assign the same spectral type but a more luminous class, B0.5Iae. 
Estimates of mi's mass lie in the range 21.7 M (van Paradijs et 
al. 1977a, assuming /=78.8°) to 24.0 M Q (Rawls et al. 2011). 
Van Kerkwijk et al. (1995) derive a value for the stellar radius, 
Ri =29.9-30.2 R Q , under the assumption that the supergiant star 
fills its effective Roche lobe at periastron. Rawls et al. (2011) 
derived/? i=3 1.82, under the same assumption. These values are 
consistent with the radii of B0.5Ia (26-38 R Q ) stars listed in 
the catalogue of Passinetti-Fracassini et al. (2001), but not with 
those of B0.5Ib (18-26 R Q ) stars. Hence, GP Vel's radius is quite 
uncertain but most likely lies in the range 26-32 R . Its effective 
temperature, r e //=25000±1000 K was determined by Sadakane 
et al. (1985) from the equivalent width of UV Fe IV and Fe III 
lines, and a comparison with other B0-B1 lb stars. Fraser et al. 
(2010) determined log(g/cm s~ 2 )=2.90 and r e , // =26500 K from 
a model atmosphere fit to the spectrum. 

The orbital inclination angle depends on GP Vel's assumed 
radius. Given the large range of possible values for the radius, it 
is, hence, rather uncertain. It also depends on the observed du- 
ration of the X-ray eclipse, 9 C , which, however, is time variable 
and energy dependent (Quaintrell et al. 2003). With these con- 
siderations in mind, most authors agree that i >78°. A recent 
analysis by Rawls et al. (2011) yields as most probable values 
/=85.9° and !=78.8°0 

Table 2 summarizes the ranges in the parameters that are 
generally adopted for the GP Vel system. The most reliable of 
its parameters are those that have been determined from the X- 
ray pulsar's delay times; specifically, P,„b, e, the projected semi- 
major axis of the neutron star orbit, o^sin/, and argument of peri- 
astron, u> per . We adopt these as the basic known parameters and 
keep them fixed throughout this paper. 

Given the considerable uncertainties that are associated with 
the remaining parameters, we opted to construct sets of self- 
consistent parameters for conducting the numerical simulations. 
The method for deriving the self-consistent parameter sets is de- 
scribed in sections 3.1-3.3. A sample of self-consistent parame- 
ters is listed in Table 3. The TIDES code was run first for a se- 
lection of these parameters. We then constructed a second block 
of parameters with small variations of one or more of the self- 
consistent parameters and the TIDES code was run for these as 
well. 



8 This is the value given in their Table 4, although in the text the value 
quoted is 77.8°. 
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Table 2. Published stellar and orbital parameters. 



Parameter 


Value 


Reference 


P (days) 


8.964368 


Bildsten et al. 1997 


e 


0.0898 


Bildsten et al. 1997 


a x s'mi(ls) 


113.89+0.13 


Bildsten et al. 1997 


9e(°) 


30-36 


van Kerkwijk et al. 1995 


t0 per ( o ) 


332.59±0.92 


Bildstein et al. 1997 


vsini (km/s) 


116±6 
56±10 


Zuiderwyk 1995 
Fraser et al. 2010 


K op , (km/s) 


17-29.7 
21.7+1.6 


van Kerkwijk et al. 1995 

"n ■ i-i ^"i r\f \ i 

Barziv et al. 2001; 




22.6+1.5 


Quaintrell et al. 2003 


log g (cm/s 2 ) 


2.90+0.2 


Fraser et al. 2010 


Teff (°K) 


26500 


Fraser et al. 2010 


25000+1000 


Sadakane et al. 1985 


V„ (km/s) 


1105 


Howarth etal. 1997 




78.8 


Rawls etal. 2011 




70.1-90. 


Quaintrell et al. 2003 


M (M \ 


24.00 


van Kerkwijk et al. 1995 
Rawls etal. 2011 




23.8 + ^ 


Barziv et al. 2001 




23.1-27.9 


Quaintrell et al. 2003 


Ri (R B ) 


31.82 


Rawls et al. 




29.9-30.2 
30.4 + '« 


van Kerkwijk et al. 1995 
Barziv et al. 2001 




26.8-32.1 


Quaintrell et al. 2003 



For the range in m\ that was analyzed, the above expression 
leads to a range R\ =27.7-28.9 R Q . It is important to note that 
R\ is the equilibrium radius of m\ in the absence of mi and as- 
suming no rotational deformation. These two effects lead to an 
aspherical shape and a unique value of the stellar radius cannot 
be defined. 



3.3. RV curve semi-amplitude, K\ 

The orbital velocity semi-amplitude is given by Binnendijk 
(1960, p. 151): 



Ki = 



27r a\ sin i 
Pil-e 2 ) 1 ' 2 



(7) 



Writing a\ = a/(l + and again using Kepler's Third Law, 



v n„r\W — -1/3 n- 1/3 (i , m i 
(1 - e 2 ) 1 ' 2 2 \ ni2 



-2/3 



(8) 



Determinations of K\ since 1976 lie in the range ~ 17-28 km/s 
(see Table 1 of Quaintrell et al. 2003 for a summary). More strin- 
gent limits have been given by Quaintrell et al. (22.6 +1.5 km/s) 
and Barziv et al. (21.7 ±1.6 km/s). These limits are derived from 
fits of Keplerian RV curves to the observational data. 



3.1. Stellar masses 



3.4. Equatorial rotation velocity 



The mass ratio, q = m2/m 1 , is constrained by the known value 
of ax sin; and Kepler's Third Law P 2 = G(m |+ m2) fl 3 ■ Using q = 
OT2//W1 = 01/02, and and a — a\ + ai — a-i(\ + q), with a\ and 
«2 = a x the semimajor axes of the two stars' orbits, 



P I Grtii sin 3 i\^ 2 ^ 
^ 2tt\ (ax sini) 3 



(4) 



Adopting P=8.964368, a x sini=H3.S9 light-seconds (Bildsten et 
al. 1997), and writing m\ in solar units, 



q = 0.22505(«« /) 3/2 (mJMo) 1 ' 2 - 1 



(5) 



Thus, for a fixed value of i, each value of m\ defines a unique 
value of m.2- The range in m\ =22.5-24.5 M leads to a range 
m 2 =1. 427-2.04 M Q for the cases analyzed for ;=78.8°. 

3.2. Primary Radius, /?i 

The value of the B-supergiant's radius, R\, is particularly rel- 
evant since the tidal forces scale as R~ 3 . There are three con- 
straints on its value. The first relies on the assumption that it 
is close to filling its Roche Lobe, an assumption that requires 
knowledge of its rotation velocity and the neutron star's mass, 
m ns . The second relies on: 1) the duration of the X-ray eclipse, 
2) the orbital separation, a, and 3) the orbital inclination, i. The 
third constraint on R\ relies on a stellar atmosphere model fit to 
its spectrum which yields log g, from which ^1 can be derived 
given knowledge of m\. We have chosen the latter constraint to 
fix R\ values for the calculations. With the observational con- 
straint on log g, the value of R t follows from the choice of m\ 

using R\ = (^y 1 ) , where g is the value of the surface gravity. 

With log g =2.90 cm s~ 2 from Fraser et al. (2010), 



There are two widely different determinations of the projected 
equatorial rotation velocity, v sin i. Zuiderwijk (1995) obtained 
v sin i=l 16 +6 km/s by fitting synthetic profiles to the observa- 
tions. This is consistent with the 1 14 km/s obtained by Howarth 
et al. (1997). Recently, however, Fraser et al. (2010) applied the 
Fourier transform method (Gray 2008) to the spectrum and ob- 
tained v sin i=56 km/s. They attribute the additional broadening 
observed in the spectral lines to "macro turbulence". 

It is important to note that, in general, the stellar rotation ve- 
locity plays a very important role in determining the behavior of 
the stellar surface. For example, faster rotation leads to a larger 
deformation of the star. In addition, the ratio of stellar angular ro- 
tation velocity, w, to orbital angular velocity, O, plays a critical 
role in the time dependence and amplitude of the tidal forcing. 
Thus, the synchronicity parameter /? = tj/Cl, which describes 
the degree of departure from an equilibrium configuration, is of 
importance in the calculation. Note that when f3-l, the system is 
in synchronous rotation, which is one of the conditions for equi- 
librium. In an eccentric binary, Q. changes with orbital phase and 
thus, /3 is a function of orbital phase. 

In our model, the stellar rotation velocity is specified as an 
input parameter through (3 Q - co/D.0, where £2 is the orbital 
angular velocity at periastron. It can be conveniently expressed 



as, 



/3o = 0.02- 



lkms-\ (1 
(P or b/days) 



ef' 2 



Ri/R 



(1 + e y/ 2 



(9) 



The two different values of v sin i imply fio ~ 0.3 or 0.6, (for 
Ri -28-29 R and i >78°). We shall refer to the case /3 Q -0.6 as 
the "standard" case and the /?o ~0.3 as the "slow" rotation case. 



fli = 5.8396 (mj/Mo) 172 ^ 



(6) 
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Table 3. Sets of self-consistent input parameters. 



ml 


m2 


i 


Ri 


v roT 


v> 

v rot 








23.6 


1.468 


78.8 


28.36 


57.09 


118.25 


0.300 


0.622 


17.37 


24.0 


1.708 


78.8 


28.60 


57.09 


118.25 


0.298 


0.617 


19.87 


24.2 


1.830 


78.8 


28.72 


57.09 


118.25 


0.296 


0.614 


21.12 


22.5 


1.427 


85.9 


27.69 


56.14 


116.30 


0.302 


0.626 


17.71 


22.7 


1.546 


85.9 


27.82 


56.14 


116.30 


0.301 


0.624 


19.02 



This table lists a sample of self-consistent input parameters derived as 
described in Section 3. Cols. 1 and 2 list the stellar masses, in M G ; 
col. 3 the orbital inclination; col. 4 the stellar radius, in Rq, cols. 5 
and 6 the "slow" and the "fast" rotation velocities, in km/s; cols. 6 and 
7 the corresponding asynchronicity parameters; and col. 8 the semi- 
amplitude of the Keplerian radial velocity curve, in km/s. 
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Fig. 1. Top: Shape of the stellar surface in the equilibrium config- 
uration (case 065); black/white color coding corresponds to the 
minimum/maximum radius; the map is centered on 180° lon- 
gitude. The asymmetry in the tidal bulges stems from the fact 
that the orbital separation is small compared to the radius of nt\ . 
Bottom: The radius as a function of longitude at the equator. 
All the 40 orbital phases for which it was computed are plotted, 
showing that no surface perturbations are present in this calcula- 
tion. 
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Table 4. Model input parameters and RV curve amplitudes. 



Case 


mi 


m 2 


i 


A 


dR/Ri 


Ri 


V 


K/,-, 




61a 


23.6 


1.468 


78.8 


0.300 


0.06 


28.363 


0.12 


18.37 


17.37 


61b 


23.6 


1.468 


78.8 


0.300 


0.10 


28.363 


0.12 


20.14 


17.37 


61c 


23.6 


1.468 


78.8 


0.622 


0.06 


28.363 


0.12 


17.70 


17.37 


61d 


23.6 


1.468 


78.8 


0.622 


0.10 


28.363 


0.22 


17.26 


17.37 


64d 


24.0 


1.708 


78.8 


0.617 


0.10 


28.602 


0.22 


20.01 


19.87 


65a 


24.2 


1.830 


78.8 


0.296 


0.06 


28.721 


0.22 


22.05 


21.12 


65b 


24.2 


1.830 


78.8 


0.296 


0.10 


28.721 


0.22 


24.37 


21.12 


65c 


24.2 


1.830 


78.8 


0.614 


0.06 


28.721 


0.22 


21.56 


21.12 


65d 


24.2 


1.830 


78.8 


0.614 


0.10 


28.721 


0.22 


21.96 


21.12 


65e 


24.2 


1.830 


78.8 


0.614 


0.08 


28.721 


0.22 


21.51 


21.12 


65f 


24.2 


1.830 


78.8 


0.614 


0.06 


28.721 


0.22 


21.34 


21.12 


71a 


22.5 


1.427 


85.9 


0.302 


0.06 


27.694 


0.12 


18.62 


17.71 


71b 


22.5 


1.427 


85.9 


0.302 


0.10 


27.694 


0.12 


20.76 


17.71 


71c 


22.5 


1.427 


85.9 


0.626 


0.06 


27.694 


0.12 


17.88 


17.71 


71d 


22.5 


1.427 


85.9 


0.626 


0.10 


27.694 


0.12 


17.99 


17.71 


72a 


22.7 


1.546 


85.9 


0.301 


0.06 


27.817 


0.22 


20.17 


19.03 


72b 


22.7 


1.546 


85.9 


0.301 


0.10 


27.817 


0.22 


22.11 


19.03 


72c 


22.7 


1.546 


85.9 


0.624 


0.06 


27.817 


0.22 


19.34 


19.03 


72d 


22.7 


1.546 


85.9 


0.624 


0.10 


27.817 


0.22 


19.38 


19.03 


72e 


22.7 


1.546 


85.9 


0.296 


0.10 


28.3 


0.22 


22.45 


19.03 


72f 


22.7 


1.546 


85.9 


0.613 


0.10 


28.3 


0.22 


19.49 


19.03 


30 


23.5 


1.44 


78.8 


0.301 


0.06 


28.3 


0.22 


18.29 


17.2 


31 


23.5 


1.44 


78.8 


0.301 


0.10 


28.3 


0.22 


20.33 


22.1 


32 


23.5 


1.44 


78.8 


0.622 


0.06 


28.3 


0.22 


17.32 


17.2 


33 


23.5 


1.44 


78.8 


0.622 


0.10 


28.3 


0.22 


17.06 


17.2 


34 


23.5 


1.44 


78.8 


0.300 


0.06 


28.8 


0.22 


18.05 


18.05 


34b 


23.5 


1.44 


78.8 


0.300 


0.10 


28.8 


0.22 


19.90 


18.05 


34c 


23.5 


1.44 


78.8 


0.620 


0.06 


28.8 


0.22 


17.55 


18.05 


34d 


23.5 


1.44 


78.8 


0.620 


0.10 


28.8 


0.22 


17.38 


18.05 


35 


24.0 


1.74 


78.8 


0.298 


0.06 


28.6 


0.22 


21.10 


20.4 


36 


24.0 


1.74 


78.8 


0.298 


0.10 


28.6 


0.22 


22.98 


20.4 


37 


24.0 


1.74 


78.8 


0.615 


0.06 


28.6 


0.22 


20.64 


20.4 


38 


24.0 


1.74 


78.8 


0.615 


0.10 


28.6 


0.22 


20.55 


20.4 


41 


24.5 


2.04 


78.8 


0.609 


0.10 


28.9 


0.22 


23.50 


23.4 


5U 


22.5 


1.52 


on n 
90.0 


U.30/ 


U.U6 


T7 h 
11 .1 


U. 12 


19.99 


19.5 


51 


22.5 


1.52 


90.0 


0.307 


0.10 


27.7 


0.12 


22.25 


19.5 


52 


22.5 


1.52 


90.0 


0.624 


0.06 


27.7 


0.12 


19.14 


19.5 


53 


22.5 


1.52 


90.0 


0.624 


0.10 


27.7 


0.12 


19.35 


19.5 


55 


23.5 


1.44 


85.9 


0.301 


0.06 


28.3 


0.12 


18.38 


18.05 


56 


23.5 


1.44 


85.9 


0.301 


0.10 


28.3 


0.12 


20.60 


18.05 


57 


23.5 


1.44 


85.9 


0.622 


0.06 


28.3 


0.12 


17.65 


18.05 


58 


23.5 


1.44 


85.9 


0.622 


0.10 


28.3 


0.12 


17.29 


18.05 



This table lists the input parameters for the models that were computed. Cases 61a-72d were computed with "self-consistent" sets of input 
parameters, as described in Section 3. m\, m 2 are given in in Solar masses; Ri in Solar Radii; v in R^/day; and K/„ and K Kep in km/s. 
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4. Results 

Table 4 lists the input parameters of the model runs performed 
for this paper. The first block of models (cases 61a-72d) are run 
with self-consistent sets of parameters, derived as is described 
in the previous section. The parameters of the second block (list 
starting with case 72e) are not fully self-consistent. Columns 2- 

5, and column 7 list, respectively, m.\ , m%, i, fio and Ri ; column 6 
contains the depth of the surface layer, dRjR\, used in the model. 
Columnn 8 lists the value of the kinematical viscosity, v (in units 
offl /day). 

4. 1 . The equilibrium case 

It is illustrative to first analyze an equilibrium case having pa- 
rameters similar to those of the GP Vel system, but in which e-0 
and (3-1. We chose the case with mi=24.2 M , m2=1.83 M , 
R\ =28.721 R and v,. 0/ =118 km/s, the latter corresponding to 
the "standard" rotation speed and an orbital inclination /=78.8°. 
In order for this system to be in equilibrium (i.e., /3=1), we had 
to set the orbital period to 12.2 days. We prefered modifying P 
instead of v ra , because in order to make f3o=l using GP Vel's or- 
bital period, we would need to set v rot ~16Q km/s, significantly 
larger than its actual rotation speed. This leads to a stronger de- 
formation of the star. The main difference of using a larger or- 
bital period instead of the actual period is that the tidal force is 
weaker, so the tidal deformation is smaller. 

The start of the TIDES calculation is characterized by a tran- 
sitory phase during which the star adjusts from its initially spher- 
ical shape to the new equilibrium shape. During this transitory 
phase, the star undergoes initially large amplitude (~0.03 R in 
the present case) oscillations that rapidly damp out, reaching an 
amplitude <0.002 R by 20 cycles after the start of the calcu- 
lation. At this time, the star has attained its equilibrium shape, 
which consists of two tidal bulges. Figure 1 shows a color-coded 
Mollweide representation of the stellar radius at each surface el- 
ement. The left edge of the map corresponds to azimuth angle 
0° (i.e., the sub-binary longitude, defined as the longitude that 
intersects the line connecting the centers of the two stars) and 
the right edge is at 360°. The plot below the map illustrates the 
magnitude of the radius at the equatorial latitude as a function 
of azimuth. Note that the size of the two tidal bulges is not the 
same due to the fact that the orbital separation is relatively small 
compared to the stellar radius. The height of the primary bulge 
above the equilibrium radius in this calculation is 6R-0.06 R . 
This corresponds to 5R/Ri =0.002, which is significantly smaller 
than the depth of the surface layer used in the TIDES calculation, 
dR/R\=0.06. Also, note that the primary bulge points directly 
towards the companion, another indication that the system is 
in the equilibrium configuration. In non-equilibrium configura- 
tions, the bulge either "lags" behind or is advanced with respect 
to the line connecting the centers of the stars, a phenomenon that 
leads to the tidal torques in the system. 

Once the equilibrium configuration is attained, there are no 
significant surface motions. Figure 1 (bottom) is actually a plot 
of the equatorial latitude of the star at the 40 orbital phases for 
which it was computed. The width of the curve is an indication 
of the stability over time of the configuration. As a consequence, 
no significant variability in the photospheric absorption-line pro- 
files is observed in our calculations^ 



9 Note that our model does not take into account a variable temper- 
ature structure across the observable stellar disk. In the case of a more 
complete calculation in which a model stellar atmosphere is used, a 
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Fig. 2. Top: Shape of the stellar surface for "standard" cases 
(model 65c) at periastron. Color coding and orientation are the 
same as in Figure 1. Bottom: The crosses represent the radius 
at each azimuth angle along the equator. The dotted line is the 
equilibrium shape shown in Figure 1, rescaled by a factor of 3. 
The strong departures from the equilibrim configuration when 
/6o=0.6 are evident. The detailed pattern changes as a function of 
orbital phase, causing the line-profile variations. 

4.2. The "standard" case 

The situation is very different when the system departs from the 
equilibrium configuration. The fact that /3o #1 implies that the 
star rotates asynchronously, and this introduces perturbations 
on the stellar surface. We illustrate this case with a model for 
mi=24.2 M , ;7i2=1.83 M , i?i=28.721 R (case 65c). Contrary 
to the equilibrium case, the stellar surface does not attain a sim- 
ple shape with two tidal bulges. The map in Figure 2 shows 
a color-coded representation of the stellar radius at each sur- 
face element, with white/black indicating maximum/minimum 
extent. The map corresponds to the time of periastron passage. 
The Mollweide representation and orientation are the same as 
in Figure 1 . The primary tidal bulge can be seen to lie between 
320° and 360°, "lagging" behind the line connecting the cen- 
ters of the two stars (which lies at azumuth 0°). This is as ex- 
pected for a sub-synchronously rotating star. Furthermore, the 
bulge shape is not smooth, but consists of smaller-scale struc- 
tures. These can be most easily visualized in the plot shown be- 
low the map, which shows the radius at the equator plotted as a 



phase-dependent variation in the line-profiles is expected due to gravity 
darkening effects, see Zuiderwijket al. 1977; and Palate & Rauw, 201 1. 



7 



Koenigsberger, Moreno & Harrington: RV curve of Vela X-l 




Radial Velocity (km s~') Phase 



Fig. 3. Predicted line profiles at different orbital phases for 
case 65c. The dotted profiles correspond to the non-perturbed, 
rigidly-rotating surface. The asymmetries in the perturbed line- 
profiles are responsible for the departures from a Keplerian RV 
curve. 



function of azimuth. A similar configuration obtains for all other 
orbital phases as well. The comparison between Figure 2 and 
the equilibrium configuration shown in Figure 1 shows that an 
"equilibrium tide" representation for GP Vel is a poor approxi- 
mation. 

The non-equilibrium configuration leads to large-scale flows 
on the stellar surface (see Tassoul 1987 and Eggleton et al. 1998 
for a good description of this phenomenon). These are referred 
to as "tidal flows". The tidal flows are motions of localized re- 
gions of the stellar surface relative to the underlying, (assumed) 
rigidly-rotating stellar interior. The changing flow patterns lead 
to variable line-profiles. A sample of the type of variability pre- 
dicted by the TIDES code is illustrated in Figure 3, where the 
perturbed line profiles are compared to their corresponding non- 
perturbed line profiles. Note the appearance of "bumps" and 
"wiggles" in the profiles, as well as the occassional blue or 
red extended wing. In calculations performed with significantly 
smaller "turbulent" speeds (for example, 10 km/s instead of the 
30 km/s used here), the "bumps" are narrower, and the profiles 
give the appearance of having narrow discrete absorption fea- 
tures that generally travel from "blue" to "red" along the line 
profile (see Harrington et al. 2009 for an example of this type 
of behavior). In the case shown in Figure 3 the discrete absorp- 
tions are significantly smoothed out due to the large "turbulent" 
speed. 

The RV curve derived from these line profiles, using the flux- 
weighted mean method, is plotted in Figure 4 (top)(3- Also plot- 
ted in this figure is the RV curve corresponding to the actual 
orbital motion (dots), and the best-fit Keplerian curve derived 
from the analysis that is described in Section 2.3 (dashes). The 
deviations of the perturbed RV curve with respect to that of the 
actual orbital motion are shown below the RV curve. They are 
as large as +5 km/s in the phase intervals </>=0. 1-0.3 and 0=0.9- 



10 A similar curve is derived if we plot the RV's obtained by fitting a 
Gaussian to the line profiles 
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Fig. 4. Radial velocity curves from the tidally perturbed line pro- 
files (squares) compared to the actual projected orbital motion 
(dotted curve) and best-fit curve (dash curve). ARV is the dif- 
ference between the perturbed and non-perturbed RV values. 
Top: m„ s =1.830 M , j6o=0.614 (the "standard" case; case65c); 
Bottom: ra M =1.468 M©, /Jo— 0.300 (case 61a). 



1 .0. Despite the large deformations on the perturbed RV curve, 
the best-fit Keplerian is very similar to the curve which describes 
the actual orbital motion. Thus, although the peak-to-peak am- 
plitude of the perturbed RV curve is significantly larger than the 
actual orbital motion, the "excursions" of the RV curve are such 
that they nearly cancel out in the fitting algorithm. This indicates 
that Keplerian curve fits to the tidally-perturbed RV curves yield 
reasonable results as long as the entire set of data is fit; i.e., with- 
out excluding portions of the RV curve. 

The general characteristics described above appear in all of 
the model runs listed in Table 4 for which /3q ~0.6. The semi- 
amplitudes of the corresponding RV curves derived from the fits 
(Ki S q) are plotted as filled symbols in Figure 5 with m.2=m ns on 
the abscissa. Also plotted in this figure are the values of K\ for 
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plitudes and, thus, it is crucial to have a well-established value 
for this parameter in order to properly estimate the magnitude of 
the RV curve perturbations. 

4.4. Dependence on layer depth 

Calculations performed with different layer depths in the "stan- 
dard" rotation case yield similar results. However, significant 
differences are present in the "slow" rotation cases where the 
dR/Ri=0.10 layer depths lead to larger RV curve amplitudes 
than those for dR/Ri=0.06. Figure 6 illustrates the radial ve- 
locities obtained from calculations with these two layer depths. 
Also plotted are the corresponding best-fit Keplerian curves and 
the curve that describes the actual orbital motion (continuous 
curve). 

The largest amplitudes occur for ySo ~0.3 and dR/R[=0.10, 
as illustrated in Figure 7 where we plot K/ sq as a function of m„ s . 
These model calculations indicate that the /3q ~0.3 systems with 
m„ s as low as 1.47 M Q are consistent with the observations. 



Fig. 5. Predicted semi-amplitude of the GP Vel RV curves for 
different assumed neutron star masses and model input parame- 
ters: open symbols: fto ~0.3; filled symbols: ySo ~0.6. All mod- 
els were run with dR/Ri-0.06. The dotted lines indicate the 
observed range of semi-amplitudes from Barziv et al. (2001) 
and Quaintrell et al. (2003). The Keplerian semi-amplitudes 
(crosses) are generally smaller than those predicted by the model 
calculations. Note that for the slow-rotation (J3o ~0.3) models, 
masses as low as m„ s =1.55 M Q yield K\ values within the ob- 
served range. 

the Keplerian orbit (crosses). The observational limits are drawn 
with horizontal dotted lines. Only models with m ns -mi >1.7 
M have values of Ki sq that lie within the observational limits. 
Hence, we conclude that for the fio ~0.6 group of models, the 
systematic shifts in RV caused by tidal flows produce only small 
departures from the actual orbital RV curve, and the derived neu- 
tron star mass is not severely affected by the tidal flows. 

4.3. The "slow" rotation case 

The slow rotation models also show significant line-profile vari- 
ability. The deformation of the RV curves leads to systemati- 
cally larger values of Ki sq . To illustrate this case we use a model 
with mi =23.6 M , m 2 =1.468 M , and Ri =28.363 Rq (case 61a). 
Figure 4 (bottom) shows that the perturbed RV curve departs 
significantly more from its Keplerian motion than the fto ~0.6 
case. There are ~5 km/s excesses around 0=0.1 and 0.65, phases 
which lie close to the extrema in the actual orbital RV curve. 
This causes the best-fit RV curve (dashes) to have a significantly 
larger amplitude than the curve describing the orbital motion 
(dots). 

The summary of Keplerian fit semi-amplitudes, Ki sq , for the 
Po ~0.3 cases is shown in Figure 5 with open symbols. For mod- 
els in which the layer depth is dR/R\-Q.Q6, we see that masses 
as low as ra Jls =1.55 yield values of K/ sq (marginally) consistent 
with the observations. 

Thus, for the ft ~0.3 group of models, the systematic shifts 
in RV caused by tidal flows produce significant deviations from 
the actual orbital RV curve. 

The basic conclusion of this section is that different rotation 
speeds lead to tidal flows with different characteristics and am- 



4.5. The "blue dip" and a structured wind 

A prominent feature in our model RV curves is a "blue dip" that 
appears at (f> ~0.3, shortly after the maximum in the curve. It is a 
persistent feature in all our computations, and it is caused by the 
asymmetrical shape of the line profiles. Figure 8 displays the line 
profiles over the phase interval 0.12-0.36 for model 65c, illus- 
trating the source of this feature. Between phases 0.14-0.22, the 
line develops an excess blue absorption wing while at the same 
time the red wing becomes less extended. The combination of 
these two effects moves the centroid of the line towards shorter 
wavelengths. A natural question is what causes this asymmetry 
in the line profiles? 

A map of the azimuthal velocity perturbations over the vis- 
ible surface of the B-supergiant at orbital phase 0=0.2 is pre- 
sented in Figure 9. The light color indicates motion of the sur- 
face elements that is faster than the underlying rigid-body rota- 
tion rate. Hence, the large white area on the left side of the map 
indicates that the surface elements that lie near the limb of the 
star are approaching the observer faster than the stellar rotation 
velocity. This leads to the more extended blue wing in the line 
profiles. Also evident in Figure 9 is the dark area in the map, 
which corresponds to surface elements whose azimuthal veloc- 
ity is slower than that of the stellar rotation. When this region 
reaches the right limb of the visible disk, it leads to a less ex- 
tended red wing in the absorption line profile. 

An extensive analysis of GP Vel's observational RV curve 
has been made by Barziv et al. (2001) and Quaintrell et al. 
(2003), using independent data sets. Barziv et al. (2001) report 
the presence of a "blue excursion" in the H<5 data that produces 
a local minimum in the RV curve at <fi ~0.37. They interpret the 
"blue excursion" in terms of a photoionization wake. Figure 10 
is a plot of the H<S 4102 A line RV values (crosses) from Barziv et 
al. (2001) showing that the "blue excursion" initiates at <p ~0.25. 
This coincides with the "blue dip" in our models, but its duration 
is far longer than that of the "blue dip". An alternative explana- 
tion for the "blue excursion" is the presence of enhanced mass 
outflow after periastron passage and it is tempting to speculate 
that it could be triggered by the strong tidal effects that appear 
after periastron passage. 

It is also interesting to note that Kreykenbohm et al. (2008) 
detected flaring activity and temporary quasi-periodic oscilla- 
tions in INTEGRAL X-ray observations. In particular, the two 
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Fig. 7. Same as in Figure 5, but from calculations performed with 
a layer depth dR/R^OAO. 
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Fig. 6. Radial velocities obtained from calculations with dif- 
ferent layer depths for the "standard" rotation velocity (top; 
cases 61c,d) and the "slow" velocity (bottom; cases 61a,b). Also 
plotted are the corresponding best-fit Keplerian curves and the 
curve that describes the actual orbital motion (continuous curve). 
Crosses and the dotted curve correspond to dR/Ri=0.06; and tri- 
angles and the dash curve correspond to dR/Ri-OAQ. The per- 
turbations due to tidal flows are more significant for the "slow" 
rotation velocity cases. 
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Fig. 8. Montage of line profiles at orbital phases </>=0. 12-0.36, 
showing that the "dip" in the RV curve is largely a consequence 
of the more extended blue wing and less extended red wing. 
Dotted curves are the unperturbed line profiles. Phases with re- 
spect to periastron passage are listed. 



largest flares observed in these data occur around orbital phase 
0.4 (with respect to periastron) in two different orbital cycles. 
This is significant because this is just ~0.03 in phase later than 
the minimum in the "blue excursion" in Barziv et al.'s data. 
Assuming a mean wind velocity of 1105 km/s (Howarth et al. 
1997), a gas stream that originates near the B-supergiant surface 
would travel ~3xl0 12 cm during this time interval, which is on 
the order of magnitude of the distance to the neutron star. Hence, 
the picture which emerges is one in which tidal flows may cause 
instabilities that produce time-dependent ouflows leading to a 
structured wind. When high-density wind regions reach the neu- 



tron star, the larger accretion rates lead to an increase in the X- 
ray emission, as described by Kreykenbohm et al. (2008). These 
authors also found temporary quasi-periodic oscillations with a 
timescale of ~2 hr which they attribute to alternating high and 
low density regions within GP Vel's stellar wind. A similarly 
structured wind has been suggested to exist in the Be/X-ray bi- 
nary system 2S01 14+650 and to possibly be produced by tidally 
induced pulsations (Koenigsberger et al. 2006). 
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Fig. 9. Map depicting the azimuthal ("horizontal") velocity field 
at orbital phase 0=0.2 for case65c. Velocities that are larger than 
the rigid-body rotation rate are coded in light color, while those 
that are slower are coded in dark. The "blue dip" in the the- 
oretical RV curve is caused by the extended blue wing in the 
photospheric absorption lines, and this is produced by the large 
white area near the approaching (left) limb of the star. The bot- 
tom panel is a plot of the azimuthal velocity perturbations along 
the equatorial latitude. 

5. Discussion and Conclusions 

In this paper, we explored the manner in which surface motions 
that are induced by the tidal interactions affect the shape of the 
radial velocity curve. The objective was to analyze the manner 
in which these motions affect the neutron star mass determina- 
tion. We performed a calculation from first principles which in- 
volves solving the equations of motions of a Lagrangian grid of 
surface elements to determine the surface velocity field on the B- 
supergiant. The velocities on the visible surface of the star were 
then projected along the line-of-sight to the observer and the lo- 
cal line-profile at each surface element was Doppler-shifted cor- 
respondingly. The Doppler-shifted local line-profiles were com- 
bined to obtain the absorption-line profile that would be detected 
by an observer at each orbital phase. For each set of model bi- 
nary parameters, the centroids of the resulting line-profiles were 
measured and a theoretical RV curve constructed. 
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Fig. 10. Theoretical RV curves for a neutron-star mass 
m ns =1.83 M Q for y8o=0.614 (continuous curve; case 65c) and 
m ns =1.468 M for/3o=0.300 (dash curve; case 61b). Data points 
are the RV measurements from Barziv et al. (2001): the average 
of several lines (triangles) and the values for H5 (crosses). The 
vertical dash lines indicate the phase of X-ray eclipse, the ap- 
proximate location of the blue "dip" in the RV curve and the 
"blue excursion". The latter may be due to outflows from the B- 
supergiant. We suggest these outflows may be associated with 
activity induced by the tidal flows. 



In all our models, the peak-to-peak amplitude of the model 
RV curve is larger than the corresponding Keplerian RV curve. 
However, the fit using Keplerian curves yields semi-amplitudes 
that may have only small departures from the true orbital motion, 
depending on the assumed stellar rotation rate. In the "standard" 
rotation case, the deformations of the RV curve are such that they 
nearly cancel out in the Keplerian fit. Thus, the derived semi- 
amplitude of the RV curve in this case is only slightly larger 
than that of the actual orbital motion. If GP Vel rotates at the 
"standard" rotation rate, then the neutron star mass may indeed 
be >1.7 M . However, in the "slow" rotation velocity cases, the 
deformations of the RV curve are significantly more pronounced 
and the Keplerian fits to the RV curves yield significanlty larger 
amplitudes than that of the orbital motion. In this case, m„ s ~ 1 .5 
Mq is feasible. 

Thus, the assumed rotation speed clearly plays a very im- 
portant role in determining the effect of the tidal flows on the 
amplitude of the RV curve. The "standard" and "slow" speeds 
we have adopted here are the two values cited in the literature. 
However, because of the prominent line-profile variability, it is 
not clear whether either of these two values accurately repre- 
sents the rotation velocity of the B-supergiant. This is an issue 
that requires further study. 

The analysis of the RV curve is further complicated by the 
observed variations that occur from one orbit to the next. Trial 
calculations do show small differences in the predicted shape 
of the star for different orbital cycles, and the maximum extent 
of the primary tidal bulge varies from one periastron passage to 
the next. However, these do not lead to significant differences in 
the predicted line-profiles at the same orbital phase for different 
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Fig. 11. The neutron star mass plotted vs. the asynchronic- 
ity factor l-/3o for the 7 binary systems listed in Table 5. 
Synchronously-rotating optical components lie at l-y8o=0. Error 
bars are the uncertainties quoted by the authors listed in Table 
5. A trend for more massive neutron stars with larger departures 
from synchronicity is observed. Whether this indicates a system- 
atic error in the mass determination due to tidal effects is an open 
question that requires a larger sample of systems before a firm 
conclusion may be drawn. 

cycles, and thus the theoretical RV curves are virtually the same 
from one orbital cycle to the next. 

A final question to address is the following: If the distortions 
in the radial velocity curve introduced by tidal flows on the sur- 
face of GP Vel lead to an artificially large neutron star mass, 
shouldn't this same effect be present in other binary systems ? 
The answer is that, indeed, the effect should be present in binary 
systems in which the stellar rotation departs from synchronous 
rotation. In this context, it is interesting to consider the set of 
neutron star binary systems with optical counterparts which is 
listed in Table 5. These are the 7 known neutron stars with op- 
tical counterparts for which stellar parameters have been well 
determined. Five of these have OB-type spectral types, while the 
remainder have A-K spectral types. The orbital periods range 
from ~1 day to ~10 days. Figure 1 1 is a plot of the derived neu- 
tron star mass vs. (1. -y8oX which tells us whether there is a de- 
pendence of m ns on the degree of departure of the optical compo- 
nent from synchronous rotation. A possible trend for larger m ns 
with larger departures from synchronicity is apparent in Figure 
11, although more data points are clearly required before we can 
conclude that a correlation exists. However, this result is consis- 
tent with the idea that asynchronous rotation leads to significant 
line-profile variability which, in turn, distorts the RV curve, po- 
tentially leading to an error in the stellar masses. Furthermore, 
this hypothesis applies to all types of binary systems. 

In conclusion, we find that: 1) the assumption of an "equilib- 
rium tide" configuration for the B-supergiant shape is a poor ap- 
proximation; 2) the tidal effects produce orbital phase-dependent 
variations in the line profiles leading to asymmetries such as ex- 
tended blue or red wings; 3) the line-profile variability causes 
the shape of the RV curve to depart significantly from that of a 
Keplerian RV curve, artificially enhancing the semi-amplitude; 



4) the characteristics of the RV curve depend on the assumed 
rotation velocity of the B-supergiant. 

Our analysis neglects the effects of a non-uniform tempera- 
ture distribution over the stellar surface, those of non-radial pul- 
sations and the possible presence of enhanced mass-outflows af- 
ter periastron passage. These effects are expected to contribute 
even further towards deforming the line profiles. Hence, we con- 
clude that given these uncertainties and the results of our calcu- 
lations, a low-mass neutron star (m ns ~1.5 M ) cannot as yet be 
excluded for the Vela X-l system. 
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Table 5. Model input parameters and RV curve amplitudes. 
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References: 'van der Meer et al. (2007) and references therein; 4 Rawls 
et al. 2011; 5 Clark 2000; 6 Mukherjee et al. 2006; 7 van Kerkwijk et al. 
1995; 8 Barziv et al. 2001; 9 Kreykenbohm et al. 2008; 10 Podsiadlowski& 
Rappaport 2000; "Casares et al. 1998. 
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